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5.1 Introduction 

For many years space researchers realized that in some specific situations solar 
dynamic power systems can provide significant savings in life cycle costs when 
compared with conventional photovoltaic power systems with battery storage 
(Menetrey 1963; Secunde et al. 1989; Prisnjakov 1991; Prisnjakov et al. 1991). A 
standard solar dynamic power system uses a mirror to concentrate solar radiation 
onto an absorber structure. By conduction through a solid material or circulation 
of a working fluid, the absorber heat is transferred to a thermal engine (i.e. a 
turbogenerator, Stirling engine, termocouple or thermionic-emitter). Alternators 
coupled to these thermal engines may generate electrical energy. Previous practice 
proved that three different cycles can be used for thermodynamic conversion of 
solar radiation: Brayton, Rankine and Stirling (Menetrey 1963; Prisnjakov 1991; 
Prisnjakov et al. 1991). For continuos operation during dark periods the use of 
melted materials to store thermal energy is being considered. Among the 
advantages of dynamic power systems one could quote their ability to provide 
electrical energy and heat simultaneously, the fact that the power plant may be 
unified by using either solar or nuclear energy or their relative invulnerability to 
corpuscular particles and to electromagnetic radiation and the possibility of power 
control according to a given power consumption schedule (Prisnjakov 1991). 

Earlier studies reported solar dynamic systems based on thermodynamic cycles 
and thermionic converters, respectively. Two systems of the former type were de¬ 
scribed in Menetrey (1963). One of them is a 3 kW solar dynamic system of 10 % 
overall efficiency consisting of a parabolic concentrator and a turbo-generator unit 
which uses mercury as working fluid in a Rankine cycle. The other one is a 15 kW 
system for space vehicle applications. Its working fluid is rubidium while the 
overall electrical efficiency is 21.7 %. The solar - thermionic system was designed 
to provide 135 W at Mars distance from the Sun with a system weighting about 
14 kg. Other studies are using permanent gases (mainly mixtures of xenon and he¬ 
lium) as working fluid in Brayton cycles. They have output power of 3 to 5 kW, 
specific weights up to 40-45 kg/kW and overall efficiency of 31.2 % (Prisnjakov 
1991; Prisnjakov et al. 1991, 1994). Significant efforts to develop solar dynamic 
power systems (25 kW) for application to a low Earth orbit space station have 
been conducted by NASA Lewis Research Center (LeRC) as an integral part of 
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the US space station development program (McLallian et al. 1988). The culmina¬ 
tion of the LeRC programme demonstration program was a full system test in a 
thermal/vacuum tank of a Space Station Freedom type solar dynamic system at the 
2 kWe size in February 1995. Related reports can be found in Secunde et al. 
(1989) and Weingartner et al. (1994). A comparison between the various power 
systems designed for the Space Station Freedom shows the sun-to-user efficiency 
is 5.7 to 6.8 % for the PV system and 18.7 to 19.8 % for the dynamic system. If, in 
addition, the storage subsystem is taken into consideration, the solar dynamic 
system becomes even more attractive, mainly because the efficiency of the latent 
heat storage is higher compared to the Ni-Cd batteries of the PV system (for 
details see Weingartner et al. (1994)). 

A natural question arising about power systems to be used on Mars surface is: 
are properly designed dynamic systems comparable in performance with the PV - 
based power systems? Some of the above references seem to indicate an 
affirmative answer. However, they all refer to systems operating in the 
interplanetary space. A solar power system placed on Mars surface has to take into 
consideration one very important feature of the Martian weather. The period be¬ 
tween the areocentric longitude L.= 161 - 326° and variants thereof have been re¬ 
ferred to as the “dust storm season”. This period is nearly centered on perihelion, 
which is the time of maximum insolation on Mars (Martin and Zurek 1993; Hour- 
din et al. 1995). In the Northern hemisphere, the Martian atmosphere varies 
greatly from year to year during the dust storm season while during the other sea¬ 
sons, most phenomena are remarkably similar from year to year. A comprehensive 
description of Martian meteorology is provided by Zurek et al. (1992). Chapter 2 
gives some information on the subject. The dust storm periods are characterized 
by an almost vanishing direct solar radiation (see e.g. Badescu (1998a) and Chap. 
2). The main consequence is that isotropic diffuse solar radiation can not be con¬ 
centrated and flat plate solar collectors have to be used in this case. 

The structure of this chapter is as follows. Upper bounds for the efficiency of 
solar thermal power plants operating in the Martian environment are evaluated in 
section 5.2. A model of a solar power plant based on a Carnot cycle is presented in 
section 5.3 while section 5.4 refers to solar Stirling engines based on selective flat- 
plate converters. These models provide lower upper limits for the performance of 
actual thermal solar engines. Section 5.5 contains the conclusions. 


5.2 Upper Limits for Solar Energy Conversion Efficiency into 
Work 

The purpose of this section is to give information about the upper limits of the 
maximum efficiency of converting solar energy into mechanical work on the sur¬ 
face of Mars. This subject was already treated in a number of papers in case of so¬ 
lar energy conversion systems placed on Earth surface. Changing the conversion 
system on Mars does not require a different theoretical framework. However, the 
characteristics of Mars meteorology is expected to have considerable influence on 
the numerical results. 
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The upper bounds for the performance presented in this section are obtained by 
using highly idealized models. They have however the advantage of giving a per¬ 
spective for the results obtained with the more realistic models developed in the 
next sections. 


5.2.1 Exergy of Solar Radiation on Mars 


The exergy of thermal radiation was studied in many papers, starting with the 
seminal work of (Petela 1964). Controversial arguments were presented in litera¬ 
ture (see e.g. Bejan (1988) for an early review). More recently, the discussion was 
re-opened in relation with the exergy of solar radiation (Petela 2003; Candau 2003). 
Also, many authors studied the thermodynamics of converting radiation energy into 
mechanical work (for early reviews see Landsberg and Tonge (1979), Bejan 
(1988). Three different theories were proposed by Jeter (1981), by Landsberg, 
Petela and Press (Petela 1961; Landsberg and Mallison 1976; Press 1976) and by 
Spanner (1964), respectively. The three theories were unified in Bejan (1987) while 
a generalization was given in Badescu (1988). The exergy factor proposed by 
Landsberg-Petela-Press, which is very often used in literature, is: 
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where T 0 and T s are the ambient and Sun temperatures, respectively. 

This factor shows the fraction of thermal radiation energy that can be converted 
into work. Note that rj LPP depends on just two temperatures (namely the tempera¬ 
tures of the two heat reservoirs). 

Let us apply the above result in case of solar radiation on Mars. The tempera¬ 
ture of solar direct radiation is T s =5760 K while the average temperature of Mars 
surface is 7’ 0 =248 K. Equation (5.1) allows to compute the value of the Landsberg- 
Petela-Press exergy factor: r\ LPP = 0.9246. For a solar direct irradiance at ground 
level (p s =400 W/m 2 , the exergy flux is 377.04 W/nr. 


5.2.2 Accurate Upper Bounds for Solar Energy Conversion Efficiency 


The theories shortly reviewed in the previous subsection predict efficiencies of 
converting radiation energy into work too high to be of practical interest. Much 
more accurate rigorous upper bounds for the real solar energy conversion effi¬ 
ciency t] were reported in recent years (Badescu 1998b, 1998c). Some of these 
new results are listed below and applied to the special case of solar energy conver¬ 
sion on Mars surface. 

(i) In case of a smaller than hemispherical blackbody radiation source and a 
blackbody receiver coupled to a reversible engine the upper bound of tj is 
(Badescu 1998b): 
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where the geometric factor 


/ 



(5.3) 


is a function of the solid angle O subtended by the radiation source. In case of a 
spherical source of radiation (e.g. the Sun) the following simple relation applies: 

/ = sin~ 1/2 S (5.4) 


where 5 is the half-angle of the cone subtended by the source. Of course, when a 
hemispherical source of radiation is considered, S = ji/2 and //, given by Eq. (5.2) 
reduces to r] LPP given by Eq. (5.1). However, generally //, < t] LPP . One can see that 
r/i depends on just two temperatures (namely the temperatures of the two heat res¬ 
ervoirs) in a way similar to the Carnot efficiency. But it also takes account of 
the geometry of the radiation source and this additional information increases its 
accuracy. 

Let us apply Eqs. (5.2)-(5.4) to solar energy conversion on Mars. The case of 
direct solar radiation is considered. Then, S is the half-angle of the cone subtended 
by the Sun when viewed from Mars. The mean distance Sun to Mars 
is dsun-Mars= 227.94-10 6 km while Sun radius is R s =696 -1 ()’ km. Consequently, S= 
tari l (R/d Sun _ M a rs )~sm l (R/d Sun -Mars)=0-003053. Equation (5.4) yields/= 18.096. Use 
of Eq. (5.2) gives the upper bound //, = 0.08393. This is obviously a more accurate 
upper bound than those predicted in section 5.2.1. 

(ii) In case of smaller than hemispherical blackbody radiation source and a 
blackbody receiver coupled to a Chambadal-Novikov-Curzon-Ahlborn (CNCA) 
engine the upper bound of // is (Badescu 1998b) 
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This is of course a more accurate estimate for // that tj I is, because part of the entropy 
generation related to heat engine operation is taken into account. As an application, 
using Eq. (5.5) in case of solar direct radiation on Mars yields ;/, 7 = 0.02542 which is 
lower than the value derived previously for //„ as expected. 

Note that a formula corresponding to t]u> P , in case the system contains a CNCA 
engine, can be obtained from Eq. (5.5) by making f = 1. One finds 
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Equation (5.6) provides a more accurate estimate for // than t] LPP , but a less accu¬ 
rate estimate than //„. Using Eq. (5.6) in case of solar direct radiation on Mars 
yields ;/', 7 = 0.72395. This is a rather high upper bound for the real efficiency and 
this loss in accuracy is the cost of using a simpler upper limit formula. 

(iii) Repeating the procedure from Bejan (1987) but including details about some 
irreversibilities one finds relationships generalizing Jeter, Spanner and Landsberg- 
Petela-Press efficiencies, respectively (Badescu 1999). For example, the relation 
generalizing the last efficiency is: 
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where n is a parameter related to the type of thermal engine. In case of n = 1 (i.e. 
a Carnot engine) one finds of course from Eq. (5.7) the efficiency ij LPP . Let us 
consider a CNCA engine (n = 1 / 2) as part of a power system converting the en¬ 
ergy of direct solar energy on Mars. Then Eq. (5.7) yields the upper limit t] ii( = 
0.76286. This rather high upper bound is comparable in size to 

(iv) The case of gray thermal radiation is of special interest, as it is sometimes 
used as an approximation for diffuse solar radiation. Note that gray radiation is a 
particular case of diluted radiation. A general theory was first developed by 
Landsberg and co-workers (see for instance the review paper by Landsberg and 
Tonge (1979)). However, important results were derived more recently and they 
will be summarized now. An accurate upper bound for the conversion efficiency 
of solar diluted radiation into work is given by Badescu (1990) 
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Here the dimensionless factor f d can be obtained from Eqs. (10) and (14) of 
Badescu (1990): 
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where Cl d is the solid angle subtended by the source of diluted radiation source 
while x(s d ) is a function of the radiation dilution factor £ d . This function was cal¬ 
culated in Landsberg and Tonge (1979). It is such that y( 1 )= I and can be ap¬ 
proximated for small £ d by: 

x(e d ) = 0.9652 -0.2777 In e d + 0.0511*,,. (5.10) 
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The above general formulation yields elegant results in the particular case of 
solar scattered radiation. The upper bound for the efficiency of converting into 
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work the energy of multiply scattered solar radiation was derived in Badescu 
(1991) and is given by 
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where the factor f di depends on the number z of scatterings (z=l,2,3,4). It is com¬ 
puted by 
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In case the scattered radiation is not subsequently concentrated, the following rela¬ 
tions apply for the dilution factors and solid angles of /-scattered solar radiation, 
respectively: 
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where £l s is the solid angle subtended by the Sun. The half-angle 8 di of the cone 
subtended by the source of /-scattered radiation can be easily computed from: 

Q-dj = 2#(l - cos Sd,i ) ■ (5.15) 


The effective temperature T di of /-scattered solar radiation is given by: 
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Table 5.1 shows results obtained in case of /-scattered solar radiation on Mars. 
One can see that after a single scattering (i.e. z'=l) solar radiation is still strongly 
anisotropic (compare 8 d ,i = 3.061° with the half-angle d s = 0.1749° of the cone 
subtended by the Sun). After four scatterings solar radiation is isotropic and an ob¬ 
server at the ground would see a uniformly brilliant sky. The effective temperature 
of the scattered radiation decreases by increasing the number of scatterings, as ex¬ 
pected. However, the upper bound efficiency rj di predicted by Eq. (5.11) yields the 
same value, whatever the number of scatterings. This is a consequence of keeping 
the energy flux constant during the scattering process. 
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Table 5.1 Upper bounds for the conversion efficiency of /-scattered solar radiation into 
work on Mars. 


Number of 

£d,i 

G c i,i 

5d,i 

T d ,i 

i)d,i 

c 

v -'max,i 

^|d,conc,i 

scatterings, i 



(deg) 

(K) 




1 

0.00000932 

0.00897 

3.06 

1377.0 

0.084 

350.70 

0.760 

2 

0.00285100 

0.41880 

21.04 

531.2 

0.084 

7.76 

0.393 

3 

0.12890000 

2.99700 

58.47 

344.8 

0.084 

1.38 

0.130 

4 

0.72650000 

6.28300 

90.00 

316.0 

0.084 

1.00 

0.084 


Sd,j - dilution factor. and 5d,i- solid angle and half-angle of the cone subtending the 
source of radiation, respectively. T di - effective temperature. r) dj - upper bound efficiency. 
Cmax.i ■ maximum concentration ratio. r|d,conc,i - upper bound efficiency for concentrated i- 
scattered radiation. Sun temperature T s =5760 K; average temperature on Mars surface 
T 0 =248 K. 


For small /-values scattered solar radiation is generally anisotropic. Consequently, 
it can be concentrated. The maximum concentration ratio C maxa is given by: 

-t 

(5.17) 

The upper bound efficiency tj,i, conC} , for fully concentrated scattered radiation is 
given again by Eq. (5.11) with f di replaced by f' di given by: 

rilki) (/ = 1,23,4). (5.18) 

7=1 

The upper bound efficiency tjd.conc.i lies between 0.760 and 0.084 for single and 
four scatterings, respectively (Table 5.1). Generally, the effective temperature of 
scattered radiation is lower, and the upper bound efficiency is higher, on Mars 
than on Earth (compare present Table 5.1 with Table 1 of Badescu (1991)). This is 
mainly caused by the lower ambient temperature on Mars. 

(v) The upper bounds at (i) to (iii) above apply to any source of blackbody un¬ 
diluted or diluted radiation. Even more accurate (but also more complex) upper 
bounds of // were derived in case of solar energy applications (Badescu 1998c, 
2000). These upper bounds do not involve the absorber temperature but are ex¬ 
pressed just in terms of the two heat reservoir temperatures. However, they take 
into account a large number of details concerning the conversion process as for 
example, the dilution factors of solar and ambient radiation, the polarization 
degree of solar and absorber emitted radiation and the thermal and optical proper¬ 
ties of the selective absorber (convective and conductive heat loss coefficients, 
concentration ratio, effective transmittance-absorptance product and emittances). 
These details increase of course the accuracy of the efficiency formula. As an ex¬ 
ample we shall use the accurate upper bound efficiency Eq. (36) of Badescu 
(1998c). For the purposes of this chapter we shall simplify it to take account just 
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on the selective properties of the absorber (i.e. its absorptance a and emittance e). 
Then, the upper bound efficiency is given by: 
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For an usual selective flat-plate solar collector with a = 0.9 and e = 0.1, Eq. 
(5.19) yields //, = 0.3724 which is considerably higher than the upper bound effi¬ 
ciency I], = 0.08393 we previously found in case of black-body (non-selective) so¬ 
lar collectors. 

Finally, note that the results presented in (i), (ii) and (v) correspond to solar 
direct radiation at normal incidence (i.e. Sun zenith angle 8 0 = 0). In the general 
case the geometric factors in Eqs. (5.3), (5.4) and (5.20) should be multiplied by 
cos'Y),, (Landsberg and Badescu 2000). For oblique radiation incidence (()<, > 0) 
this procedure will yield lower estimates for the upper bound efficiency. 


5.3 Carnot Solar Power Plants with External Irreversibility 

In Sect. 5.2 we presented simple accurate upper bound formulas for the efficiency 
of solar power plants operating on Mars. Much more involved models are neces¬ 
sary to further increase the accuracy of results. A solar power plant model based 
on the Carnot cycle with external irreversibility is presented here. This approach 
keeps the advantage of generality but also takes into account the irreversibilities 
associated to the heat transfer at the hot and cold components of the plant. 

5.3.1 Solar Power Plant Model 

Generally, a solar power plant consists of a solar collector, a thermal engine and a 
heat storage system. A solar collector mainly consists in mirror, absorber (re¬ 
ceiver), transparent cover(s) and thermal insulation. The mirror is missing in case 
of flat - plate collectors while for some concentration systems the receiver is not 
protected by transparent cover(s). The solar power plant analysed in this section is 
designed to operate during the whole Martian year. Consequently, it consists of a 
flat-plate selective solar collector coupled to a thermal engine. No heat storage 
system is considered in this approach. 

5.3.1.1 Solar Collector Model 

A selective flat - plate solar collector is shown in Fig. 5.1. Its main components 
are: the absorber plate (1) with pipes for the working fluid and a transparent cover 
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(i) (ii) ( m ) 


Fig. 5.1 Solar collector, (i) longitudinal section; (ii) transversal section; (iii) lateral view. 
1 - absorber plat; 2 - carbon dioxide layer; 3 - transparent cover; 4 - bottom thermal insulation 


(3) of thickness a. The gap between (1) and (3) (thickness s ) is filled with carbon 
dioxide at the Martian atmospheric pressure (carbon dioxide represents 95 % of 
the Martian atmosphere). The bottom thermal insulation (4) has a thickness b. The 
collector tilt angle is [i. One denotes by T c and T a the absorber plate and ambient 
temperatures, respectively. Other temperatures related to collector operation are 
defined in Table 5.2. 

Table 5.2 Various temperatures related to solar collector operation 


Zc 

T si 


T 

± SS 

T 

J a 


T Qs ,=(Tc+T si )fl 


T s ,ave-(T s i+T ss )/2 
Tss,a = (T S s + T a )/2. 


T lz =(T c +T a )/2 


absorber plate temperature 

temperature of lower surface of transparent cover 
temperature of higher surface of transparent cover 
atmosphere temperature 

average fluid temperature between absorber plate and 
lower surface of transparent cover (see Fig. 5.2) 
average transparent cover temperature 
average temperature of the atmosphere boundary 
layer near the collector transparent cover 
average bottom insulation temperature 


Direct beam (b), diffuse (d) and ground reflected (r) fluxes of solar energy are 
incident on the collector. Their magnitude depends upon collector's tilt and orien¬ 
tation. Generally, the incident and absorbed solar energy density fluxes ( <p inc and 
<p abs , respectively) are given by 


Vine ~ Gb +Gd + Gr 

Vabs = ( Ta )b Gb + {ra) d (Gd + G ,) • 


(5.21,22) 


Here G 's denote solar irradiances at collector level while (ra)'s are effective 
transmittance - absorptance products. 

The collector heat losses are studied by using the thermal resistance method 
(Eaton and Blum 1975) (see Fig. 5.2 and Table 5.3). Then, the flux of heat losses 
from the collector towards the ambient is given by; 
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R i 


R 2 


r 4 


R 5 



Fig. 5.2 Solar collector thermal resistances. T a - ambient temperature; T c - absorber plate 
temperature; T si - temperature of lower surface of transparent cover; T ss - temperature of 
upper surface of transparent cover 


Table 5.3 Thermal resistances associated to solar collector operation (see Figs. 5.1 and 5.2) 


Thermal Process 

Resistance 

Ri Convection and conduction between absorber plate 

(1) and transparent cover (3) 

R 2 Radiation between absorber plate (1) and the lower surface of 

the transparent cover (3) 

R 3 Conduction through transparent cover (3) 

R 4 Radiation between the higher transparent cover surface (3) 

and Martian atmosphere 

R 5 Convection between the higher transparent cover surface (3) 

and Martian atmosphere 

R 6 Conduction through the bottom thermal insulation(4) 


q tot = Ul(t c - Ta)= ( Tc - To)/ R t „t 


(5.23) 


where Cl is the overall heat losses coefficient while R, ot is the total thermal resis¬ 
tance. Cl is given by (see Fig. 5.2): 
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where the equivalent resistance Rj. 5 is computed from: 
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Details concerning the evaluation of thermal resistances follow. The thermal resis¬ 
tance R 1 is computed from 
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Table 5.4 (a) Relationships for convection heat transfer between absorber plate and trans¬ 
parent cover (Eaton and Blum 1975); (b) General relationships for convection heat transfer 
between two plane parallel surfaces (Stefanescu et al. 1982, p. 160) 


(a) 

N uc,si - UH + N uv 



N UH = 0.195 -Gr 025 

for 2 • 10 3 <Gr <4 ■ 10 5 


N u H = 0-068 Gr 0 - 33 

for Gr>4-10 5 


N u v = 0.180 Gr°- 25 {slLf 11 

for 1 • 10 4 <Gr<2- 10 s 


N u v = 0.065 Gr 0J3 {s/Lf 11 

for 2- 10 5 <Gr <11- W 6 

(b) 

N uc,si-l 

for Gr ■ Pr<l ■ 10 3 


N uc ,si = 0.105 {Gr Prf 30 

for 1 • 10 3 <Gr ■ Pr<l • 10 6 


Nu c , si = 0.400 (Gr-Prf 20 

for 1- 10 4 <Gr<2-10 5 


N uc ,si = 0.180 {Gr Prf 25 

for Gr ■ Pr> 1 ■ IQ 10 


H and V denote a horizontal and vertical solar collector; /3 [rad] and L are tilt angle and 
collector length, respectively. 


where <lco 2 >C,si is thermal conductivity of carbon dioxide and Nuc, S i is Nusselt 
number, both of them evaluated at the average fluid temperature T Csi . 

To compute Nuc.si, a number of criterial relationships developed for solar col¬ 
lectors are used (Table 5.4a) (Eaton and Blum 1975). For cases not covered by 
Table 5.4a we use the general relationships from Stefanescu et al. (1982) (Table 
5.4b). Prandtl and Grashoff numbers in Table 5.4 were computed with: 

p r _ < c P.C02 > C.si < V C o 2 >C,si _ Sm PcOi S ^ C ~ T si) (5 ? 7 ) 

^ ^-CO-2 >C,si < 0 ("a , K( (),TC,si 

where g M is gravitational acceleration on Mars while Pco 2 and R COl are the pres¬ 
sure and gas constant for the carbon dioxide between absorber plate and transpar¬ 
ent cover. In Eq. (5.27) the specific heat at constant pressure (< c p COl > C si ), the 
dynamic viscosity (< > C si ) and the thermal conductivity (< A COl > Csi ) of 
carbon dioxide were interpolated at the average temperature T Csi . In computations 
the thermal properties of carbon dioxide from Pop et al. (1987) were used. 

The thermal resistance R 3 is given by: 


R3 = a/A s - (5.29) 

The temperature dependence of transparent cover's thermal conductivity (/',) was 
described by using an interpolation parabolic function whose coefficients were ob¬ 
tained by using experimental data (the transparent cover is made of crystal; thermal 
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properties are found in Lide (1991)). The average transparent cover temperature 
T save is used in computations. 

The thermal resistance R 4 is evaluated with: 



where T ssa is the average temperature of the atmosphere boundary layer near the 
collector transparent cover. 

Because of the wind, the heat transfer between the transparent cover and Mar¬ 
tian atmosphere is mainly by forced convection. Thus, the thermal resistance R s 
was computed from: 



where L is collector's length on wind direction and < A COl is carbon dioxide 
thermal conductivity at temperature T ssa . In Eq. (5.31) Nu ssa depends on the criti¬ 
cal length X T which separates the laminar and turbulent wind flows over the col¬ 
lector. The relationships of Table 5.5 are used. The Reynolds and Prandtl numbers 
are evaluated with: 


Re _ Pcq 2 WL p r _ < C P.C0 2 >™-g K ’ho 2 >ss - a 5 32 ) 

< 7 h:o 2 >ss,a Rco 2 T ss,a < A COi ^SSyd 

where W is wind speed and the thermal properties of the atmosphere (< ij (:0l > ss a , 
< c P .co-, >ss,a and< A COl > ss a ) are interpolated at the average temperature T ssa . 

The thermal resistance R 6 is evaluated with: 

R 6 = b/A iz (5.33) 

where is bottom insulation thermal conductivity, computed at the average tem¬ 
perature T iz . 


Table 5.5 Relationships for forced convection heat transfer between transparent cover and 
atmosphere (Stefanescu et al. 1982, p. 126) 


„ _ <r 7 „5 <r, C02 >ss - a R C0 2 Tssa 

X-T - 2 ' 10 il) 

_ Pco, w _ 

N u S sa = 0.335 ■ Re 0 ' 50 ■ p r 0 50 
N u SS4 , = 0.036- Re 0 ' 80 ■ p r 075 


for L< Xt (laminar flow) 
for L> Xt (turbulent flow) 
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5.3.1.2 Engine Model 

A variety of working fluids were proposed to be used in space solar engines. They 
include mercury, potassium and rubidium in the early studies (Menetrey 1963) or 
organic fluids and mixtures of noble gases more recently (Prisnjakov et al. 1991; 
Angelino and Invernizzi 1993). The later seem to be more appropriate in case of 
Martian solar power plants. The largest and smallest working fluid temperatures 
are denoted T' and T" , respectively (Fig. 5.3). 



Fig. 5.3 Temperature - entropy diagram for a Carnot cycle with external irreversibility (for 
notations see text) 

Generally, T c > T' and T" > T„ ■ Hence, both 

X = T C -T' and Y = T" - T a (5.34) 

are positive quantities. The flux of solar energy (f) h absorbed by the collector is 

<t>abs = (pabsAc , (5.35) 

where A c is collector surface area. Part of this flux is transferred to the working 
fluid ( Qc ) while the other part constitutes the flux of thermal losses to the ambi¬ 
ent ( Qi ). Then, the solar collector steady-state energy balance is: 

</> a b s -Q c -Q L = Q • (5-36) 

The heat fluxes Q c and Q L can be written as: 

Q c =hcAcX and Q L = U L Ac(Tc~T a ) > (5.37) 

where h c and U/ are appropriately defined overall heat transfer coefficients be¬ 
tween the solar collector and the working fluid and ambient, respectively. 

The Carnot engine partially converts the flux Qc into power. Thus: 
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where W is the output power and Q„ is the thermal energy flux finally reaching 
the surroundings. There is no further increase of entropy during work production. 
Consequently 

Q£ + Ql = o (5.39) 

T T" 

The heat flux Q a can be written as 

Q a =h a A a Y , (5.40) 

where A a and h a are appropriately defined engine radiator heat transfer area and 
overall heat transfer coefficient between working fluid and ambient, respectively. 


5.3.1.3 Performance Indicators 


Four performance indicators may be considered in relation to Martian solar plant 
operation (Badescu et al. 1999, 2000a, 2000b). First, there is the Carnot engine ef¬ 
ficiency defined as: 

T1 . =1-^—. (5.41) 

• engine v ' 

Second, there is the output power W , evaluated from: 

W=h c A c X-h a A a Y. (5.42) 

Flere we used the Eqs. (5.37), (5.38) and (5.40). Third, there is the system effi¬ 
ciency, defined as: 

W 

system ~ ~ ' (5.43) 

(PabsAc 

Finally, the solar (or sun-to-user) efficiency is defined as: 


W 

$ absA-C 


(5.44) 


Information about some of these indicators is reported in the following. The exer- 
getic efficiency has been not reported for reasons of space. 


5.3.1.4 Optim isation 

The output power W is maximised now. The weight of space power plants should 
be kept at minimum. Following previous studies (see e.g. Badescu (1994)) we ex¬ 
press the mass of the total heat transfer area (i.e. the mass M c of the solar collector 
and the mass M a of the engine radiator) as a A: fraction from the total mass M lol of 
the space system: 
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Mc+M a —kM, ot . (5.45) 

One denotes d c and d a the superficial mass density of collector and radiator, re¬ 
spectively. One can simply write: 

Mc—dcAc > Ma — daAa- (5.46) 

One can use dc=d a as a reasonable assumption. By using Eqs. (5.45) and (5.46) 
one derives the following constraint to be fulfilled: 

A c + A a = A = k^£L. (5.47) 

d a 

Other constraints (say F i and F 2 ) are the Eqs. (5.36) and (5.39). By using Eqs. 
(5.34), (5.37) and (5.40) they become: 

P ~ hc Ac X hg (A - Ac) _ q 
' Tc' x Ta + Y (5.48,49) 

F 2 = (p abs -U L {Tc-Ta)-h c X = 0 

The power W is considered here as a function of T c , X, Y and A c . The La¬ 
grange function L associated to W is: 

L = W + A,Fi + A 2 F 2 , (5.50) 

where li and X 2 are multipliers. The maximum output power is obtained by solv¬ 
ing the six equations shown in Table 5.6. 

Table 5.6 Equations to be solved for power maximisation (for details see text) 


1. 

2 . 

3. 

4. 

5. 


-=-a / , cAc * -^°o 

(Tc-X f 

= he Ac + Ai j C - A 2 he = 0 

(Tc-Xf 

--hAA-Ae)-x M A ' Ac) J‘' =0 

( T a + Y f 

-=h c X+ ha Y+A -^7 + -^] = 0 

VTc-X Ta-YJ 
_ he Ac X h a ( A - Ac )Y _ „ 


Tc-X 


Tg + Y 


-=(Pabs-UL{Tc-T a )-h c X = 0 


6 . 
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Table 5.7 Optimum parameters for maximum output power (for details see text) 


1. 


3. 


Tc, 


opt ' 


1 + f- 


1+f- U+f) 


K 


1/ 2 


IT 


+ K 

1/ 2 


1 / 2 


Ul 


X opt 


_(l + f )-(l + f )" ‘ K t 
he ° 


Y opt = 


27 2 ( 7 + ^ 


4. Ac.opi 


1 + K 

K 1/2 

1 + K 


1/2 


T a 


1/2 


Table 5.7 shows the solution that applies for K > 1. There, the dimensionless 
parameters K and / are defined as: 


K= — 
he 


K ( H h c 


1 + - 


V 


U 


f = 


LJ 


(Pabs 

YhJa 


(5.51) 


Table 5.7 and Eq. (5.42) allow writing the maximum power as: 

• h a A 2 

W mm =—Y 0 ~ pt (5.52) 

The optimum solar efficiency r] S0 i ari0pl can be derived by using Eqs. (5.44), (5.52) 
and Table 5.7. 


5.3.2 Results and Discussions 

One of the goals of this chapter is to perform a rough comparison between the per¬ 
formance of PV cell power systems and dynamic solar power plants. In case of PV 
cell systems we had in mind the Pathfinder's Sojourner. Sojourner is a small (11.5 
kg), six - wheel robotic vehicle built at Jet Propulsion Laboratories. Sojourner was 
landed on Mars aboard the Pathfinder spacecraft on July 4, 1997. In the same sol 
(i.e. Martian day) she began to traverse the Martian terrain, perform science and 
technology experiments, and transmit images and data back to the Lander space 
craft. Sojourner's equipment (computers, motors, radio modem) was mainly pow¬ 
ered by a lightweight 0.34 kg solar array of receiving surface area 0.22 m 2 . The 
PV cell system was designed to provide Sojourner with around 16 W of electric 
power at noon on Mars. The landing place was in Ares Vallis at 19.17° N and 
33.21° W on the U.S. Geological Survey (USGS) cartographic network (Golom- 
bek et al. 1997). Pathfinder landed in late northern summer (areocentric longitude 
L v =143°) and operated for 83 sols. 

To allow comparison with Sojourner performance an appropriately designed 
dynamic power plant equipped with a flat - plate selective solar collector is ana¬ 
lyzed here (see Table 5.8). A number of assumptions are accepted as follows. A 
simple relation is used for the effective transmittance - absorptance products: 
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Table 5.8 Details about selective flat - plate solar collector design 


Quantity 

Notation 

Value 

Number of transparent covers 

N 

1 

Transparent cover thickness 

a 

0.003 m 

Distance between transparent cover and absorber plate 

s 

0.045 m 

Bottom thermal insulation thickness 

b 

0.1 m 

Short wavelengths transparent cover transmittance 

Tcs 

0.82 

Short wavelength transparent cover absorptance 

a ps 

0.90 

Long wavelengths (IR) absorber plate emittance 

Sl 

0.10 

Long wavelengths (IR) transparent cover emittance 

£2 

0.88 

Thermal conductivity of bottom insulation (polyure¬ 
thane) 

Ab 

0.02 W/(mK) 


(A c =0.22 m 2 ) 


(t a) b = (t a) d = T cs a ps , (5.53) 

where z cs and a ps are the transparent cover transmittance and absorber absorptance, 
respectively, both of them for short (solar) wavelengths. The collector's transpar¬ 
ent cover is made up of crystal (light flint glass). Its thermal conductivity lies be¬ 
tween 0.691 W(mK)" 1 at -100 °C and 1.025 W(mK)' 1 at +100 °C (Lide 1991). The 
ratio 0.34/11.5=0.0296 between the PV array weight and total Sojourner weight 
can be used as a first estimate for the coefficient k appearing in Eq. (5.47). Pre¬ 
liminary tests were performed to determine an optimum collector area A Copt close 
to Sojourner PV cell surface area (0.22 m 2 ). Finally, the value A =0.23 m 2 was 
used in computations. It corresponds to a radiator superficial mass density d a = 
1.48 kg/m 2 in Eq. (5.47). It is rather close to the value d a = 2.3 kg/m 2 accepted in 
Angelo and Buden (1991) in case of the radiator of a nuclear power satellite but 
only half of the value d B =3 kg/m 2 accepted in Mozjorine et al. (1991) for a solar 
space power station. 

The distance between VL1 and Pathfinder Lander is around 815 km. Thus, one 
expects quite similar meteorological and actinometric features in both places. 

Four strategies of collecting solar energy were considered in preliminary tests 
(Badescu 1998a). The following two strategies will be used in this section: (i) 
horizontal collector - strategy H; (ii) the collector tilt and orientation are continu¬ 
ously adjusted to keep the receiving surface perpendicular on Sun's rays - strategy 
P. Strategy H is easier to use while strategy P gives the higher power output most 
time of the year. 

5.3.2.1 Solar Collector Operation 

Solar collector thermal losses mainly depend on the heat transfer regime between 
absorber plate and transparent cover (separated through the layer 2 of thickness s 
in Fig. 5.1). We considered the case of a horizontal solar collector. Three design 
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Table 5.9 Three design solutions for the system solar collector - thermal engine 


Case 

ha 

[ W/(m 2 K) ] 

he 

[ W/(m 2 K) ] 

Thermal engine performance 
level 

(a) 

1 

1 

Low efficiency 

(b) 

10 

10 

Medium efficiency 

(c) 

500 

50 

High efficiency 


h a = heat transfer coefficient between the working fluid passing through the radiator of the 
heat engine and ambient (see Eq. (5.40)); h c = heat transfer coefficient between solar collec¬ 
tor and working fluid (see Eq. (5.37)). 

solutions for the system solar collector - thermal engine were analysed (Table 5.9). 
They correspond to three different performance levels. For analysis we have cho¬ 
sen the middle of a clear sky autumn sol, when the thermal losses are rather high 
(sol 301 VL1, 14.04 local solar time (LST)). [We remind that sol numbering 
started when each Viking Lander touched down the Martian soil.] 

The dependence of the overall heat loss coefficient U L on the distance s be¬ 
tween absorber plate and transparent cover is shown in Fig. 5.4 for the three cases 
of Table 5.9. The heat losses decrease by increasing s . This is specific to the con¬ 
duction heat transfer regime. The thickness s is limited, however, by economic 
reasons. In the following we use ,?=0.045 m and the high efficiency engine (case 
(c) in Table 5.9). Note that the chosen value of.? is about two times larger than the 
usual value for Earth placed solar collectors. 

For Earth placed solar collectors the heat transfer between absorber plate and 
transparent cover is mainly by convection. This prompted various technical solu¬ 
tions meant to keep convection at local level (the honeycombs and Francia cells 
are examples of solutions for convection suppression) (see Meinel and Meinel 



Fig. 5.4 Overall heat loss coefficient t/ L as a function of parameter s (distance between ab¬ 
sorber plate and transparent cover - see Fig. 5.1). Three design solutions (a), (b),(c) are con¬ 
sidered (see Table 5.9). (VL1 site, autumn, sol 301, 14.04 LST) 
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(1976, p. 404)). The convective heat transfer is activated when the product 
Gr -Pr exceeds 1000 (see Table 5.4b). We analysed the heat transfer above the 
absorber plate by using the whole set of available meteorological data (Table 2.2 
of Chap. 2). The results show that the product Gr ■ Pr varies between 9.65 and 
211.2, well below the threshold value of 1000. We draw the important conclusion 
that on Mars the thermal losses between absorber plate and transparent cover are 
exclusively by conduction. 

The convection heat losses from transparent cover to atmosphere are mainly 
controlled by wind speed. The change from laminar to turbulent convection occurs 
in case the heat transfer surface exceeds a certain critical length X f (see Table 5.5). 
The flow regime above the transparent cover was analysed by using the whole set 
of available meteorological data. Results are shown in Fig. 5.5. The minimum 
value of the critical length X T is 5.43 m. Consequently, the convection heat trans¬ 
fer above the transparent cover is laminar. The critical length X r is smaller during 
winter and larger during summer (Badescu et al. 2000a). No significant difference 
between VL1 and VL2 sites is observed (note that a larger number of meteoro¬ 
logical values are available for VL2 - see Table 2.2 of Chap. 2). The critical length 
obviously decreases when the optical depth increases. During the "clear sky" days 
(T < 1 ) the critical length values have an important dispersion. Scattering dimin¬ 
ishes during the dust storm days (T> 1) . 

For the design solution we selected (case (c) in Table 5.9), the values of the over¬ 
all heat loss coefficient JJ L range between about 0.5 and 1 Wm" 2 K _1 (Fig. 5.6). 
Broadly, these values are comparable with those of Earth located vacuum solar 



Fig. 5.5 Critical length X T (see Table 5.5) as a function of optical depth. A high efficiency 
engine was considered (see Table 5.9). All the available meteorological data were used. 
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Fig. 5.6 Overall heat loss coefficient t/ L as a function of optical depth. A high efficiency 
engine was considered (see Table 5.9). All the available meteorological data were used. 


collectors (see e.g. Meinel and Meinel (1976, p. 387)). The dispersion of the L\ val¬ 
ues is slightly smaller during winter and higher during autumn (Badescu et al. 
2000a). No obvious dependence of U\ on the optical depth is observed. However, 
the dispersion of the U L values is larger for values r < /, i.e. during the "clear" days. 
Also, Fig. 5.6 shows that U L does not depend significantly on the latitude of the so¬ 
lar collector. 

5.3.2.2 Thermal Engine Operation 

The parameters /z a and h c have an important influence on performance. For analy¬ 
sis we have chosen the middle of a clear sky autumn day (sol 301 VL1, 14.04 
LST). Some results are shown in Fig. 5.7 for the maximum power provided by the 
thermal engine. The influence of both h. d and h c is more important on the maxi¬ 
mum power provided by the heat engine that on the overall efficiency (Badescu 
et al. 2000a). 

Generally, the influence on performance is strong for low values of both heat 
transfer coefficients /z a and h c (less than 100 WnrK" 1 . The performances are 
practically not dependent on the higher values of h c . The influence of h d on per¬ 
formance is more significant. The strategy of collecting solar radiation has a rather 
weak influence on the dependence of performance on h. d and h c . As expected, 
strategy P leads to higher performance (Badescu et al. 2000a). 

The maximum output power could be as high as 12 W and 13 W in case of 
strategies H and P of collecting solar radiation, respectively (Badescu et al. 2000a). 
This is lower than the designed power of Sojourner PV arrays (16 W at noon). 
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Fig. 5.7 Isolines of maximum power Wmax provided by the thermal engine as a function of 
the overall heat transfer coefficients h d and h c . Horizontal solar collector (strategy H). Un¬ 
concentrated solar radiation was considered (VL1 site, autumn, sol 301, 14.04 LST). 


However, this comparison put in disadvantage the solar dynamic system as we used 
input meteorological data from late autumn at VL1 site while Sojourner was de¬ 
signed to operate during Martian summer. This is a consequence of the fact that 
none of the 22 complete records available during summer at VL1 site (see Table 
2.2 of Chap. 2) corresponds to solar noon. 

5.3.2.3 Performance Dependence on Atmospheric Optical Depth 

Previous results showed that the dynamic solar power plant performance depends 
significantly on the heat transfer coefficients h. d and h c . Practical situations should 
normally lye between the cases (b) (i.e. medium efficiency engine) and (c) (i.e. 
high efficiency engine) in Table 5.9. These two cases will be used to provide 
lower and upper bounds for actual performance. To test the performance of the 



Optical depth 

Fig. 5.8 Dependence of solar efficiency ^ so iar, 0 pt on the atmospheric optical depth at VL2 
site. High efficiency engine (see Table 5.9) and horizontal solar collector. All the available 
meteorological data for VL2 site were used. 
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thermal engine - solar collector combination we simulated its operation at both 
VL1 and VL2 sites. We analyzed both the horizontal collector (strategy H) and the 
collector oriented perpendicularly on sun's rays (strategy P), respectively (Badescu 
et al. 2001a). At VL1 site we used meteorological data from summer, autumn and 
winter, year 1. All available data were used to simulate the operation at VL2 site. 
Some results are shown in Fig. 5.8 (for solar energy conversion efficiency) and 
Fig. 5.9 (for output power). 

The influence of latitude on performance is obvious. When the collector is 
horizontal, the solar energy conversion efficiency is generally smaller at VL1 site 
as compared to VL2 site (Fig. 5.8). When a P collector is considered, the effi¬ 
ciency increases (at small T 's) but it still is smaller at VL1 site that at VL2 site 
(Badescu et al. 2001a). The maximum power provided by a horizontal collector at 
VL1 site shows some interesting features (Badescu et al. 2001a). Generally, it is 
smaller than in case of VL2 site. However, there are some situations (for T ~ I) 
when the meteorological effects compensate the latitudinal effects and the output 
power is quite similar for both VL1 and VL2. If a P collector is used, the power 
increases for small T 's but slighter at VL1 site than at VL2 site (Fig. 5.9). 



0 12 3 4 


Optical depth 

Fig. 5.9 The same as Fig. 5.8 in case of the maximum output power Wmax • High efficiency 
engine and P solar collector. 

In case the medium efficiency engine is coupled to a horizontal collector, the 
efficiency does not exceed 0.13 at VL2 site (Badescu et al. 2001a). Generally, it is 
lower during summer, when the horizontal orientation is surely a good option 
(Badescu 1998a). It is higher during winter, when the optical depth is higher, indi¬ 
cating dust storms occurrence. During dust storms the incident solar radiation is 
mainly diffuse and the horizontal orientation is the best strategy for radiation col¬ 
lection (Badescu 1998a). In case a medium efficiency engine is coupled to a P col¬ 
lector, the maximum efficiency values lye around 0.13 (Badescu et al. 2001a). The 
efficiency increases at small optical depth values, which generally correspond to 
summer. The centre of the output data cloud is placed around 0.09-0.10. The P 
strategy does not diminish the performance at larger optical depth. 
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For a horizontal collector attached to a high efficiency engine, the solar effi¬ 
ciency is as high as 0.18 (Fig. 5.8). The qualitative features pointed out previously 
maintain. The lower sun-to-user efficiency corresponds to lower optical depth. If a 
P collector is considered, the solar efficiency increases during summer and spring 
but does not significantly exceed 0.18 (Badescu et al. 2001a). The centre of the 
output data cloud is placed around 0.15, higher than in case of the horizontal col¬ 
lector. Again, the P strategy does not diminish the performance. This strategy 
keeps constant the solar efficiency at higher optical depths but improves the per¬ 
formance as compared to the H strategy at smaller optical depths. 

The power provided by a system consisting in a horizontal collector and a me¬ 
dium efficiency engine does not exceed 7 W (Badescu et al. 2001a). Generally, it 
is lower for small optical depth. Indeed, the flux of direct solar energy is higher at 
smaller optical depth and the horizontal orientation is not recommended in this 
case. At higher optical depth direct solar radiation diminishes and diffuse radiation 
increases. As a result, the horizontal orientation is close to optimum. This explains 
the higher values of W max around x = 1. For further increase of x the diffuse radia¬ 
tion flux diminishes. Consequently, W max decreases in this case. If a P collector is 
considered, the solar efficiency is higher during summer and spring, when the op¬ 
tical depth is small (Badescu et al. 2001a). Generally, the efficiency increases as 
compared to the H strategy. At high optical depth the solar efficiency is compara¬ 
ble for both H and P strategies. 

Using a better engine with a horizontal collector leads to output power up to 13 
W for x around 1 (during sols of spring, autumn and winter) (Badescu et al. 
2001a). For small values of x the output power is smaller. The P strategy is rec¬ 
ommended in combination with high efficiency engines (Fig. 5.9). In this case the 
sun-to-user efficiency during sols with small optical depth increases strongly, ex¬ 
ceeding 25 W. The centre of the output data cloud is placed around 13 W. The P 
strategy can be used at high optical depth, too. 

53.2.4 Performance Dependence on the Level of Solar lrradiance 

The dependence of solar efficiency on the level of incident global irradiance was 
studied in (Badescu et al. 2001a). In case of a horizontal collector operating at 
VL2 site, a slightly non-linear relationship exists between solar efficiency and in¬ 
cident global irradiance for both the medium and the high efficiency engines. 
When the P collector is considered, the dependence is more complicated. 

When a combination medium efficiency engine - horizontal collector is consid¬ 
ered, the efficiency is higher (whatever the incident irradiance is) during autumn 
and winter and lower during summer. This rather surprising fact is valid for the P 
collector, too. It is probably due to the influence of the ambient temperature, 
which is lower during autumn and winter. The situation is different in case of the 
high efficiency engine. Then, the solar efficiency is higher during summer, even if 
it corresponds to lower incident irradiance. This remark is even more obvious 
when a P collector is considered. A more uniform solar efficiency values distribu¬ 
tion is obtained in this case. 
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In case of a horizontal collector, the efficiency increases by increasing the ir- 
radiance level, whatever the type of engine is. For the P collector, the efficiency 
depends on the irradiance level in a more interesting manner. Thus, the lower 
bound of the efficiency increases with the irradiance level, but the upper bound 
practically does not depend on the irradiance level. This is valid for both the me¬ 
dium and high efficiency engines. 

During autumn and winter the dependence of ;/ s0]ar on the level of incident 
global irradiance is nearly similar for both the horizontal and the P collector oper¬ 
ating at VL1 site. Important differences exist during summer. The P collector is 
more effective, as expected. The influence of irradiance level on solar efficiency 
decreases drastically beyond a certain threshold value (about 300 W/m 2 ). This ap¬ 
plies to both types of solar collectors. 

Figure 5.10 shows some results for the output power at VL1 site. Meteorologi¬ 
cal data from year 1 were used. In case of a horizontal collector the dependence of 
the maximum output power on the level of global irradiance is not differentiated 
upon season (Badescu et al. 2001a). When a P collector is considered, the same 
incident irradiance value leads to an output power obviously higher during sum¬ 
mer. For a horizontal collector the output power increases when the input irradi¬ 
ance increases (Badescu et al. 2001a). This remark maintains for the P collector 
except for summer results (Fig. 5.10). 



Solar global irradiance [W/m 2 ] 


Fig. 5.10 Dependence of maximum output power W ma x on global irradiance at VL1 site. 
High efficiency engine (see Table 5.9) and P solar collector. Meteorological data for sum¬ 
mer, autumn and winter year 1 were used. 

The power provided by a system based on a horizontal collector increases at 
VL2 site with the level of solar irradiance, as expected (Badescu et al. 2001a). The 
increase is almost proportional in case of the medium efficiency engine and 
slightly non-linear for the high efficiency engine. Generally, there is no obvious 
dependence on season. In case of the P collector, the output power obviously de¬ 
pends on season, with a maximum during spring and summer and a minimum dur¬ 
ing autumn and winter. This is true for both the medium and high efficiency en¬ 
gines (Fig. 5.11). 
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Fig. 5.11 Dependence of the maximum output power Wmax on the level of solar global ir¬ 
radiance incident on a horizontal surface at VL2 site. High efficiency engine (see Table 5.9) 
and P solar collector. All the available meteorological data for VL2 site were used. 


The seasonal influence is similar for both solar energy conversion efficiency 
and output power in case of a medium efficiency engine - horizontal collector 
combination (Badescu et al. 2001a). This is still valid when the collector is con¬ 
nected to a high efficiency engine. However, in this case the power values have a 
larger dispersion, which increases slightly with increasing the irradiance. 

In case of a medium efficiency engine and a P collector, the influence of the 
season is stronger on power than on solar efficiency (Badescu et al. 2001a). This is 
even more obvious if one connects a P collector to a high efficiency engine. The 
high solar efficiency values during spring and autumn at smaller irradiance values 
are not the main cause of the higher output power (which is obtained in summer 
(Fig. 5.11)). For a medium efficiency engine, a P collector is recommended 
mainly during summer and spring (which are dust storm - free seasons) (Badescu 
et al. 2001a). One finds higher output power values as compared to the case when 
a horizontal collector is used. If the summer and spring values would be neglected, 
one can see that choosing a P instead of an H collector leads to a relatively small 
increase in power during the other two seasons (autumn and winter). This is valid 
for the high efficiency engine, too (Fig. 5.11). However, the increase in power by 
using a P collector during seasons without dust storms is more spectacular. Com¬ 
pare the power values as high as 20 W obtained at lower input irradiance by using 
a P collector (Fig. 5.11) with the power values of at most 12-13 W obtained at 
high irradiance values by using a horizontal collector (Badescu et al. 2001a). 

53.2.5 Diurnal Variation of Performances 

The input data file contains only a few sols with complete records covering appro¬ 
priately the whole daylight time. The most part of the sols has one to three records 
only associated to time periods close to sunrise, noon or sunset. In particular, no 
summer sol with complete records covering the whole day was found either for 
VL1 or for the VL2 site. For each season, we selected those sols that are richer in 
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available data. In case of VL1 these sols are 301 (autumn) and 328 (winter). Note 
that no available data refers to spring (see Table 2.2 of Chap. 2). In case of VL2 
the sols selected are 420 (spring) and 872 (autumn). The sol 406 (winter) is also 
available but no reference to it will be made here, as the winter conditions are rela¬ 
tively similar for both Viking Lander sites. 

Various efficiencies defined in Section 5.1.3 were reported in Badescu et al. 
(2000b) in case of sol 301 (Martian autumn) at VL1 site. During most of the time 
the optimum engine efficiency f/ eng i ne , op t lies around 33 %. The efficiency is de¬ 
creasing suddenly near sunset. Generally, the optimum engine efficiency is less 
dependent on the type of engine and on the strategy of collecting solar radiation. 
In exchange, the system efficiency ii sys tem,opt is strongly dependent on engine qual¬ 
ity and less dependent on the strategy of solar collection. It lies around 0.2 for the 
high efficiency engine and around 0.15 for the medium efficiency engine. When 
the high efficiency engine is considered, the decrease in performance is more 
abrupt near sunset. 

The optimum solar efficiency >/ so iar,opt is smaller than the optimum system effi¬ 
ciency (Fig. 5.12) but exhibits the same qualitative features. Generally, it doesn't 
exceed 0.15 for the high efficiency engine and 0.10 for the medium efficiency en¬ 
gine. The strategy of collecting solar radiation has not a strong influence on 
f/soiar.opt- In exchange, the strategy significantly influences the maximum output 
power Wmax (Fig. 5.13). If the optimum solar efficiency is almost constant during 
the day, the maximum power’s time variation is closely related to the temporal 
variation of solar irradiance. This is more obvious for the high efficiency engine. 
The maximum output power can exceed 15 W for the high efficiency engine and 
doesn't exceed 8 W for the medium efficiency engine. This proves that properly 
designed thermal power plants are comparable in performance with PV-based 
power systems (one reminds that Sojourner’s design output power is 16 W at solar 
noon). Moreover, the previous analysis puts in disadvantage the solar dynamic 



Local solar time 

Fig. 5.12 Optimum solar efficiency f7 so iar,opt during sol 301 (autumn) at VL1 site. Medium 
efficiency engine (ME) and high efficiency engine (HE) (see Table 5.9). Strategy H - hori¬ 
zontal collector; Strategy P - collector permanently kept perpendicular on Sun’s rays. 
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Fig. 5.13 The maximum output power Wmax during sol 301 (autumn) at VL1 site. For de¬ 
tails see Fig. 5.12. 

system as we used here input meteorological data from late autumn at VL1 site 
while Sojourner was designed to operate during Martian summer. 

The hourly variation of solar collector’s optimum temperature is shown in Fig. 
5.14. It is less dependent on the strategy of collecting solar radiation and a little 
more dependent on engine quality. The high efficiency engine leads to a smaller 
optimum collector temperature. Generally, the collector temperature exceeds with 
more than 100 K the ambient temperature. The hourly variation of 7c i0p t is 
stronger than that of the ambient temperature. This proves that T c , 0 p t is mainly 
controlled by the level of incident solar radiation. 

The optimum temperature of the working fluid in the cold part of the engine 
(7” opt ) is well correlated with the time variation of the ambient temperature (Badescu 
et al. 2000b). 7” opt is smaller in case of the high efficiency engine, when the influ¬ 
ence of the strategy of collecting solar radiation is obvious. This influence is less 



Fig. 5.14 Optimum solar collector temperature T c , op t associated to solar thermal power plant 
operation during sol 301 (autumn) at VL1 site. For details concerning the power plant see 
Fig. 5.12. 
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important in case of the medium efficiency engine. In exchange, the 
optimum temperature of the working fluid in the worm part of the engine (7”’ opt ) is 
better correlated with the hourly variation of the solar radiation (Badescu et al. 
2000b). Both the engine quality and the strategy of collecting solar radiation are 
important factors. Generally, 7”’ opt is close to 7 Cpipl . 

Sol 328 belongs to winter at VL1 site. The main performance indicators of the 
solar thermal power plant are studied in Badescu et al. (2000b). The optimum solar 
efficiency values are rather close to the autumn values. The quality of the engine is 
very important and the strategy of collecting solar radiation is less important, as far 
as their influence on the performance indicators is considered. This is due to the 
fact that sol 328 belongs to a dust storm period (Badescu 1998a). The maximum 
output power is, however, much smaller than during autumn. At noon it reaches 10 
W in case of the high efficiency thermal engine but doesn't exceed 5 W in case of 
the medium efficiency engine. The maximum output power time variation is 
stronger than that of solar efficiency but this is less obvious than in autumn. The 
daily variation of the optimum collector temperature is rather weak. Generally, the 
high efficiency engine leads to a lower optimum collector temperature. Also, the 
optimum temperature is lower for a collector perpendicular on Sun's rays. 

Sol 420 belongs to springtime at VL2 site. In this case, the optimum solar effi¬ 
ciency depends significantly both on engine quality and on the strategy of collect¬ 
ing solar radiation (Badescu et al. 2000b). In case of a P collector, the optimum 
solar efficiency is nearly constant during the day. It varies rather much for the 
horizontal collector (with a maximum at noon). The optimum solar efficiency val¬ 
ues are quite close to those from autumn and summer at VL1 site. The maximum 
output power shows features similar to those of the optimum solar efficiency. 
Generally, the P collector assures a more constant maximum output power during 
the day. The engine quality is very important. The high efficiency engine con¬ 
nected to a P collector can ensure a maximum output power around 15 W during 
most part of the day. This is close to the design power of the PV system of Mars 
Pathfinder’s Sojourner. When a P solar collector is considered, the optimum col¬ 
lector temperature shows a remarkable constancy during the day. 

Sol 872 corresponds to autumn at VL2 site. The optimum solar efficiency 
shows an interesting variation during the day. This is more obvious for a P collec¬ 
tor (Badescu et al. 2000b). In case of the horizontal collector, the optimum solar 
efficiency has a maximum around the noon. The optimum solar efficiency is 
higher in case of the high efficiency engine, as expected. It ranges from 8 to 16 %. 
The maximum output power shows similar features. Generally, the better per¬ 
formance in the morning corresponds to the stronger solar irradiance during that 
part of day (Badescu 1998a). Note that a horizontal collector assures a constant 
power during the day. The high efficiency engine should be used in combination 
with a P collector. A maximum output power of 15 W can be provided only a 
short time period during the day. In case of a P collector, the optimum temperature 
is rather constant until 14.00 LST. When a horizontal collector is considered, the 
optimum temperature has a maximum reached between 12.00 and 14.00 LST. 
Generally, 7’ C opt is lower than in the cases above and this can be correlated to the 
smaller ambient temperature value. 
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5.3.2.6 Correlation between Power and Conversion Efficiency 

The dependence of the maximum output power on the optimum solar efficiency at 
VL1 site is shown in Badescu et al. (2000b). Meteorological data from year 1 were 
used there. Let us consider the case of the medium efficiency engine. The maxi¬ 
mum output power increases when the optimum solar efficiency increases, as ex¬ 
pected. There is little dependence on season, except for a few summer values in 
case of the P collector, when the solar efficiency is obviously higher. Broadly 
speaking, there is little difference between the maximum output power provided 
by systems based on H and P collectors, respectively. The highest optimum solar 
efficiency is 0.11 and the maximum output power doesn’t exceed 10 W. When the 
high efficiency engine is considered, the influence of optimum solar efficiency on 
the maximum output power is important. The highest optimum solar efficiency is 
0.16 and the output power can be as high as 15 W for an H collector and 17 W for 
a P solar collector. 

Figure 5.15 shows some results obtained at VL2 site (Badescu et al. 2000b). All 
the available meteorological data for the two years of VL2 operation were used in 
computations. Generally, the dependence of the maximum output power on the 
optimum solar efficiency is non-linear. In case of a horizontal collector, the de¬ 
pendence seems to indicate a quadratic dependence (Wmax 00 Hsoiar.opt)- The coef¬ 
ficient of proportionality is smaller in case of the medium efficiency engine 
(Badescu et al. 2000b). When the P collector is considered, the dependence 
W maxC^ solar,opt) is more complicate, but keeps the quadratic feature. The P collec¬ 
tor has a more compact dispersion of the optimum solar efficiency values. It 
ranges from 8 to 12 % in case of the medium efficiency engine and between 13 
and 19 % in case of the high efficiency engine. 



Conversion efficiency 

Fig. 5.15 Dependence of the maximum output power Wmax on the optimum solar effi¬ 
ciency ;/ S0 i ari0pt at VL2 site. High efficiency engine (see Table 5.9) and horizontal solar col¬ 
lector. Meteorological data for summer, autumn and winter (year 1) and spring, summer 
and autumn (year 2) were used. 
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In all cases a certain optimum solar efficiency threshold (around 5 %) must be 
exceeded in order the system provide useful power. The dependence of the maxi¬ 
mum output power on optimum solar efficiency is stronger in case of the high ef¬ 
ficiency thermal engine than in case of the medium efficiency engine, on one 
hand, and in case of the P collector than for the horizontal collector, on the other 
one hand. It seems no obvious difference exists between the performances of the 
power plant in the two years of VL2 operation. 


5.4 Solar Stirling Power Engine 

A more realistic solar thermal engine is considered in this section. It consists of a 
solar collector - Stirling engine combination. The Stirling engine is indeed very at¬ 
tractive due to its operation at low temperature difference between the two heat 
reservoirs. This could enable power production during the Martian dust storm pe¬ 
riod, when both the incoming solar energy flux and collector temperature are 
small. Preliminary results were briefly reported in Badescu et al. (2001b) where 
the diurnal solar plant operation is studied. One concluded that the efficiency of a 
usual (Earth based designed) Stirling solar engine diminishes with about 8 % 
when operates under Martian weather conditions. However, the recommended 
volumetric ratio lies between 1 and 2.6 and the best thermal agent is helium, in 
agreement with conclusions based on Earth made experiments. The model previ¬ 
ously proposed in Badescu et al. (2001b) is now presented. A number of im¬ 
provements are also included (Badescu 2004). 

5.4.1 Solar Engine Model 

The solar engine consists of a selective flat-plate solar collector coupled to a Stir¬ 
ling engine with partial heat regeneration. No heat storage system is considered in 
this approach. 

5.4.1.1 Solar Collector Model 

First, we shall consider the flat-plate solar collector (Fig. 5.16). Its main compo¬ 
nents are the collector plate and a transparent cover of thickness a . The gap of 
thickness s between these two components is filled with carbon dioxide at Mar¬ 
tian atmospheric pressure. Solar direct and diffuse radiation penetrates the trans¬ 
parent cover and is absorbed by the collector plate. The thermal energy generated 
in the solar collector plate is transferred to the thermal agent at the hot head of the 
Stirling engine and later-on part of it is converted into work. 

Steady state energy balance for the solar collector yields 

Q„ f , s -Q'-Q l( , t =0- (5-54) 

Here Q a i, s is the thermal energy flux generated into solar collector plate, Qi ost is the 
heat flux lost towards the ambient and Q' is the heat flux transferred to the 
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engine's hot head 


Fig. 5.16 The solar collector of the Stirling engine 


thermal fluid inside the Stilling engine. The thermal energy flux generated into the 
solar collector plate is given by 

Qabs = (PabsA-c , (5.55) 


where Ac is the area of collector plate surface. The heat flux lost towards the 
ambient is given by the usual Newton relationship: 

Qiost = U l A c (T c —T a ), (5.56) 


where U L is the overall heat loss coefficient, while T c and T. d are collector and 
ambient temperature, respectively. Note that T c is an averaged value over collector 
surface and normally is not constant in time. 

Seven temperatures are associated to solar collector operation (Badescu 2004). 
With a single except (i.e. the average temperature of bottom insulation) they have 
the meaning explained in Tables 5.2 and 5.3. The collector heat losses could be 
studied by using the thermal resistance method. Then, the flux of heat losses qi ost 
per unit collector surface area is defined by: 

qiost = ^f L = U L {T c -T a ) = , (5.57) 

4 - R to , 


where R tol is the total thermal resistance given by 

r„=-^ + r, + -^L. 
Ri +R 2 R 4 + R 5 


(5.58) 


Note that R tot here is equivalent to R i_ 5 in the model of solar collector of section 
5.1.1. Therefore, Fig. 5.2 and Table 5.3 show the meaning of the thermal resis¬ 
tances R t to R 5 . Details concerning the computation of these thermal resistances 
are similar to those presented in Section 5.1.1 and will be not repeated here. 
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5.4.1.2 Stirling Engine Model 

During the 1980s research was started on low temperature differential Stirling 
engines. A first Stirling engine running with a temperature difference below 20 
degrees Celsius has been developed (Senft 1996). Meanwhile a rather large 
number of research or commercial Stirling engines using unconcentrated solar 
radiation were built, with powers from a few to tens of Watt (Badescu 2004). 

A simple Stirling engine model is used at this stage of the analysis. It was first 
developed in Howell and Bannerot (1977) and sligthly improved in Badescu 
(1992). The model is based on two simplifying assumptions. First, one neglects 
the heat losses during the thermal transfer from the solar collector to the working 
fluid inside the engine. Second, one neglects the friction of the moving parts. 
These make of course the results reported below to be upper bounds for the 
performance of a real solar Stirling system. The Stirling cycle efficiency i] Sljrlmg is 
defined as: 


t] Stirling 


P 

Q' 


7/ max 

l + £>/7max 


(5.59) 


where P is the output mechanical power while // max =1-7), IT C is the maximum 
efficiency of the cycle in case of perfect heat regeneration. In case the thermal 
agent is an ideal gas the parameter D is given by Howell and Bannerot (1977), 
Badescu (1992): 

D _ x (5-60) 

(k -1) In e v 


Here x( e [0,1]) is the heat regeneration factor, k is the adiabatic exponent of the 
thermal agent and E v is the volumic ratio (i.e. the ratio of the extreme values of the 
thermal agent total volume), x = 0 means ideal heat regeneration. In this case the 
efficiency of the Stirling engine Eq. (5.59) equals the Carnot efficiency. 

Use of Eqs. (5.54), (5.55), (5.56) and (5.59) yields the following dependence of 
the power provided P on the collector plate temperature T c \ 

P(Tc ) = — tyab,-U l(T c -T c )]A c . (5.61) 

1+ U(l-la He) 


Note that U\ is not a constant but depends on T c (among other variables) (see Fig. 
5.17). Consequently, a simple numerical procedure to maximize directly the 
power P given by Eq. (5.61) was used here, as follows. One monotonously 
increases the plate temperature T c starting from the ambient temperature. For any 
value of T c one evaluates first the overall heat loss coefficient U\ (Eq. (5.57)) and 
then the output power P (Eq. (5.61)). The maximum power F max and the 
appropriate optimum plate temperature 7c,opt are finally selected from the series of 
results. 
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Collector temperature (K) 


Fig. 5.17 Dependence of overall heat loss coefficient Ul on collector temperature. All 
available meteorological data were used (see Table 2.2 of Chap. 2). 

Figure 5.17 shows that the overall heat loss coefficient U L increases by increas¬ 
ing the collector temperature. The values of U L range between 0.3 and 0.7 Wm’ 2 K _1 . 
Broadly, these values are smaller than those of Earth located flat-plate vacuum 
solar collectors (see e.g. Benz and Beikircher (1999) where the experimentally de¬ 
rived overall heat loss coefficient is about 1.25 Wm' 2 K _1 during a rather heavy utili¬ 
zation as steam production). They are slightly smaller than the results obtained dur¬ 
ing simulation of the Carnot solar engine in section 5.3 ( Ul values between 0.5 and 
1 WnrK' 1 ). This is to be expected as an (ideal) Carnot engine normally operates at 
higher collector temperature than the Stirling engine considered here. The depend¬ 
ence of Ul on collector temperature justifies a posteriori the numerical optimization 
procedure we used in this section. 

The way of using the output power of the Stirling engine is not relevant for this 
study as the key aspect here is the comparison with the output power of a similarly 
sized PV-based system. Consequently, electric generators and/or gear boxes are 
not considered, neither efficiencies for these (possible) components are included 
in the model. 

Use of the output power P allows to define the solar energy conversion 
efficiency (or solar-to-user efficiency) as ij=P/(<p mc A c ). The maximum solar 
efficiency is given by 

Hsol, max ^max /(<PincAc) . (5.62) 

Note that in the case analysed here both Tj so i and P have their maximum for the 
same value Tc^pt- 

5.4.2 Results and Discussions 

A horizontal solar collector is considered. This is the easiest-to-use strategy. It 
gives for many periods of time a performance less than 10 % smaller than the 
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strategies involving more sophisticated, orientable, solar collectors (Badescu et al. 
2000a, 2000b, 2001a). Details about solar collector design are given in Table 5.8 
(where the information about the bottom insulation should be omitted of course). 
The values of the effective transmittance- absorptance products were computed 
with (ra) dir =(ra) d if=ra. 

Three Stirling engine thermal agents are usually considered: carbon dioxide 
(adiabatic exponent k = 1.33), air (k = 1.4) and helium (k = 1.66). Helium gives the 
best performance and is our choice for this section. A simple analysis proves that 
the Stirling engine efficiency increases by increasing the compression ratio e v - 
Here the value e v =2 is adopted as a good compromise between better perform¬ 
ance and smaller engine size. Increasing the heat regeneration parameter x from 
the ideal value x = 0 diminishes the engine output power. Here the value x = 0.2 is 
accepted, in agreement with current terrestrial practice. 

Results of computations performed by using the meteorological data during the 
autumn and winter of the first year at VL1 site are shown in Fig. 5.18. It is obvi¬ 
ous that both the maximum power P max and the maximum solar efficiency t] so i max 
increase around the solar noon. However, r/ so t max seems to be less dependent on so¬ 
lar time than P max . The maximum power of the solar Stirling engine around the 
noon is comparable in magnitude with the designed power of Sojourner’s PV cells 
system (i.e. 16 W). Also, the solar efficiency is as large as 18 %, which is close to 
the design value of Sojourner’s PV cells efficiency (about 18 %). Note, however, 
that the horizontally placed collector is a good solar energy collection strategy es¬ 
pecially during autumn and winter. Indeed, results not shown here prove that dur¬ 
ing spring and summer the solar Stirling engine performance is slightly worse. 
During the warm season other strategies of collecting solar energy are much more 
effective (Badescu 1998a). 

The large dispersion of the results in Fig. 5.18 is mainly due to the influence of 
optical depth. The performances of solar Stirling engine obviously decrease by 



Fig. 5.18 Dependence of maximum power P max and maximum solar efficiency ^ so iar,max ° n 
local solar time (in Earth hours). Results obtained by using the meteorological data during 
the autumn and winter in the first year at VL1 site. 
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increasing the atmospheric optical depth (Badescu 2004). The maximum effi¬ 
ciency yjsoimax is less dependent on the optical depth than the maximum power 
P m ax- Note than under-unitary optical depth values are usually associated with 
“clear sky” conditions while during the dust storm the optical depth has larger 
values (Badescu 2001). 

The dependence of the maximum output power P max on optimum solar effi¬ 
ciency sol,max is shown in Fig. 5.19. All the meteorological data of Table 2.2 in 
Chap. 2 were used. The maximum output power increases when the optimum solar 
efficiency increases. Generally, the dependence of P max on r] so i,„, a x seems to indi¬ 
cate a quadratic relationship. This is in agreement with previous results obtained 
in case of solar Carnot engines. 

A certain optimum solar efficiency threshold (around 5 %) must be exceeded in 
order the system provide useful power. Again, this confirms previous results re¬ 
ported in section 5.3. The optimum solar collector temperature T Copt increases by 
increasing the solar efficiency, as expected. The larger dispersion of the T c opt val¬ 
ues for the same value of the solar efficiency is mainly due to the difference in the 
solar time. 
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Conversion efficiency 


Fig. 5.19 Dependence of the maximum power P max and optimum collector temperature 
Tc.opt on maximum solar efficiency ^ so i imax . All available meteorological data were used (see 
Table 2.2 of Chap. 2). 


5.5 Conclusions 

The chapter focuses on solar thermal power generation on Mars surface. A "dy¬ 
namic" solar power plant (which consists in a solar collector - thermal engine 
combination) is proposed as an alternative for the more usual photovoltaic cells. 

Upper bounds for the efficiency of solar thermal power plants operating in the 
Martian environment are evaluated in section 5.2. The earlier theories predict too 
high efficiencies for radiation energy conversion into work to be of practical interest. 
Much more accurate simple upper bounds were already proposed in literature. For 
reader convenience they are summarized in section 5.2 and applied to power plants 
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operating on Mars under both direct and diffuse solar illumination. Both black-body 
and selective absorbers were considered. 

A more elaborated but still idealized model uses a Carnot cycle to describe so¬ 
lar engine operation in section 5.3. The solar power plant is designed to operate 
during the whole Martian year. Consequently, it contains a selective flat - plate so¬ 
lar collector. A detailed model of collector heat losses towards the ambient is 
developed. The optimization procedure developed in section 5.3.1 is based of fi¬ 
nite-time thermodynamics methods. The following two strategies of collecting so¬ 
lar energy were analysed: (i) horizontal collector - strategy H; (ii) the collector tilt 
and orientation are continuously adjusted to keep the receiving surface perpen¬ 
dicular on Sun's rays - strategy P. The first strategy is easier to implement and use 
in practice while the later strategy provides higher collected solar energy most 
time of the year. A medium and a high efficiency thermal engine were analysed. 
Practical situations should normally lye between these two cases. All the computa¬ 
tions were performed for a solar collector comparable in size with Pathfinder’s So¬ 
journer (solar energy collection area of 0.22 m 2 ). During the application we used 
as input data values of atmospheric pressure and temperature, wind speed and at¬ 
mospheric optical depth measured at Viking Landers sites. The main results are: 

(1) The heat losses between the solar collector absorber plate and transparent 
cover are exclusively by conduction. This differs from the Earth-based solar col¬ 
lector where the main heat loss mechanism is convection. 

(2) The convection heat transfer above the transparent cover is laminar. This 
differs from the Earth-based solar collectors where turbulence is the main convec¬ 
tion mechanism. 

(3) For the design solution we selected, the values of the overall heat loss coef¬ 
ficient t/ L range between about 0.5 and 1 Wm 2 K _1 . Roughly, these values are 
comparable with those of Earth located vacuum solar collectors. The dispersion of 
the U] values is slightly smaller during winter and higher during autumn. No ob¬ 
vious dependence of U\, on the optical depth is observed. However, the dispersion 
of the U] values is larger for values during the "clear" days. U\ does not depend 
significantly on the latitude of the solar collector. 

(4) Two parameters are used to quantify the heat transfer inside the thermal en¬ 
gine. They are h Ll (the heat transfer coefficient between the working fluid passing 
through the radiator of the heat engine and ambient; see Eq. (5.40)) and h c (the 
heat transfer coefficient between solar collector and working fluid (see Eq. (5.37)). 
The influence of /? a on performance is more significant. 

(5) The strategy of collecting solar radiation has a rather weak influence on the 
dependence of performances on h. d and h c . As expected, strategy P leads to higher 
performances. 

(6) The dynamic solar power plants equipped with selective flat - plate collec¬ 
tors could provide power comparable to that of similar-size PV cell systems. 

(7) The influence of latitude on performance is obvious. Generally, the solar 
efficiency is smaller at VL1 site as compared to VL2 site for both strategies of 
collecting solar radiation. In most cases the maximum power provided by a hori¬ 
zontal collector at VL1 site is smaller than at VL2 site. However, in some 
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situations the meteorological effects compensate the latitudinal effects and the 
output power is quite similar at both VL1 and VL2 sites. 

(8) In case a medium efficiency engine is coupled to a horizontal collector, the 
solar efficiency does not exceed 0.13 at VL2 site. Generally, it is lower during 
summer and higher during winter's dust storms. In case the medium efficiency en¬ 
gine is coupled to a P collector, the solar efficiency increases during summer. 

(9) The solar efficiency is as high as 0.18 in case of a horizontal collector at¬ 
tached to a high efficiency engine. If the high efficiency thermal engine is con¬ 
nected to a P collector the solar efficiency increases significantly during summer 
and spring but does not exceed 0.18. 

(10) Using a high efficiency engine coupled to a horizontal collector leads to 
output power up to 13 W during spring, autumn and winter. 

(11) The P strategy is recommended mainly during summer and spring in com¬ 
bination with high efficiency engines. In this case the output power could be as 
high as 25 W. 

(12) During an autumn day at VL1 site (sol 301) the optimum engine efficiency 
is around 33 %. The system efficiency lies around 20 % for the high efficiency 
engine and around 15 % for the medium efficiency engine. The optimum solar ef¬ 
ficiency doesn't exceed 15 % and 10 %, respectively, in the two cases. The maxi¬ 
mum output power can exceed 15 W and doesn't exceed 8 W for the high and 
medium efficiency engine, respectively. 

(13) During a winter dust-storm day at VL1 site (sol 328) the maximum output 
power is much smaller than during autumn. At noon it reaches 10 W in case of the 
high efficiency thermal engine but doesn't exceed 5 W in case of the medium effi¬ 
ciency engine. 

(14) During a spring day at VL2 site (sol 420) the optimum solar efficiency 
values are quite close to those from autumn and summer. The high efficiency en¬ 
gine connected to a solar collector kept perpendicular on Sun’s rays can ensure a 
maximum output power around 15 W during most part of the day. This is close to 
the design power of the PV system of Mars Pathfinder’s Sojourner. 

(15) During an autumn day at VL2 site (sol 872) the optimum solar efficiency 
ranges from 8 to 16 % in case of a good quality engine. The high efficiency ther¬ 
mal engine should be used in combination with a solar collector kept perpendicu¬ 
lar on Sun’s rays. A maximum output power of 15 W might be provided a short 
time period during the day. 

(16) When a horizontal solar collector is considered, the dependence of the 
maximum output power on optimum solar efficiency seems to be quadratic at both 
VL1 and VL2 sites. When a collector perpendicular on Sun’s rays is considered, 
this dependence is more complicate, but keeps the quadratic feature. 

(17) A certain optimum solar efficiency threshold (around 5 %) must be ex¬ 
ceeded in order the system provide useful power. 

(18) The dependence of the maximum output power on the optimum solar effi¬ 
ciency is stronger in case of a high efficiency thermal engine than in case of a me¬ 
dium efficiency engine, on one hand, and in case of a collector kept perpendicular 
on Sun’s rays than in case of a horizontal collector, on the other one hand. 
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(19) No obvious difference exists between power plant performances in the two 
years of VL2 operation. 

A solar Stirling engine based on a horizontal selective flat-plate converter was 
analyzed in section 5.4. A numerical optimization procedure was used to maxi¬ 
mize the power provided by the engine. The main results are as follows: 

(1) The values of the overall heat loss coefficient (J L range between about 0.3 
and 0.7 WmV. 

(2) The solar energy conversion efficiency at noon could be as high as 0.18 at VL2 
site during autumn and winter. Generally, it is lower during spring and summer. 

(3) The output power at noon produced by the Stirling engine at VL2 site was 
as high as 16 W during autumn and winter. It obviously decreases during the dust 
storm periods. 

(4) Properly designed solar Stirling engines equipped with horizontal flat-plate 
collectors could provide output power comparable in magnitude to that of PV cell 
systems. However, this conclusion should be treated with caution as it is based on 
figures referring to an existing PV system, on one hand, and to a simplified solar 
Stirling engine model, on the other one hand. 

(5) To increase the output power and the efficiency, the design parameters of 
Stirling engines operating on Mars should be different from those of terrestrial en¬ 
gines. For example, one expects the compression ratio to exceed the usual value 
e v =2 while a smaller regeneration factor than .r = 0.2 to be implemented. 
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